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Abstract 



»~ The problem of the nonequivalence of the sets of equilibrium points and energy- Casimir extremal 

points, which occurs in the noncanonical Hamiltonian formulation of equations describing ideal 

fluid and plasma dynamics, is addressed in the context of the Euler equation for an incompressible 

inviscid fluid. The problem is traced to a Casimir deficit, where Casimir elements constitute the 
00 ; 

center of the Lie-Poisson algebra underlying the Hamiltonian formulation, and this leads to a study 
of the symplectic operator defining the Poisson bracket. The kernel of the symplectic operator, for 



in 

this typical example of an infinite-dimensional Hamiltonian system for media in terms of Eulerian 



variables, is analyzed. For two-dimensional flows, a rigorously solvable system is formulated. The 
nonlinearity of the Euler equation makes the symplectic operator inhomogeneous on phase space 
(the function space of the state variable), and it is seen that this creates a singularity where 
the nullity of the symplectic operator (the "dimension" of the center) changes. Singular Casimir 
elements stemming from this singularity are unearthed using a generalization of the functional 
derivative that occurs in the Poisson bracket. 
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I. INTRODUCTION 

Equations that describe ideal fluid and plasma dynamics in terms of Eulerian variables are 
Hamiltonian in terms of noncanonical Poisson brackets, degenerate brackets in noncanonical 
coordinates. Because of degeneracy, such Poisson brackets possess Casimir elements, invari- 
ants that have been used to construct variational principles for equilibria and stability [1]. 
However, early on it was recognized that typically there are not enough Casimir elements to 
obtain all equilibria as extremal points of these variational principles. In [2, 3] it was noted 
that this Casimir deficit is attributable to rank changing of the operator that defines the 
noncanonical Poisson bracket. Thus, a mathematical study of the kernel of this operator 
is indicated [4], and this is the main purpose of the present article. In order to facilitate 
this study it is necessary to place the noncanonical Hamiltonian formalism on a more rig- 
orous footing, and this subsidiary purpose is addressed in the context of Euler's equation 
of fluid mechanics, although the ideas presented are of more general applicability than this 
particular example. 

We start by reviewing finite-dimensional canonical and noncanonical Hamiltonian me- 
chanics, in order to formulate our problem. These dynamical system have the form 

f t =Jd z H(z), (1) 

where z = (z 1 , z 2 , . . . , z m ) denotes a set of phase space coordinates, H is the Hamiltonian 
function with d z its gradient, and the mxm matrix J (variously called e.g. the cosymplectic 
form, Poisson tensor, or symplectic operator) is the essence of the Poisson bracket and 
determines important Lie algebraic properties [2, 5]. For canonical Hamiltonian systems of 
dimension m = 2n the matrix J has the form 

Jc= I (2) 

Noncanonical Hamiltonian systems allow a z-dependent J(z) (assumed here to be a holo- 
morphic function) to have a kernel, i.e. Rank(J(z)) may be less than m and may change as 
a function of z. 

From (1) it is evident that equilibrium points of the dynamics, i.e. points for which 
dz/dt = Vt, satisfy 

d z H(z) = 0. (3) 



However, in the noncanonical case these may not be the only equilibrium points of a given 
Hamiltonian H(z), because degeneracy gives rise to Casimir elements C(z), nontrivial (non- 
constant) solutions to the differential equation 

J(z)d z C(z) = 0. (4) 

Given such a C(z), replacement of the Hamiltonian by H(z) + C{z) does not change the 
dynamics. Thus, an extremal point of 

d z [H(z) + C(z)]=0 (5) 

will also give an equilibrium point. Note, in light of the homogeneity of (4) an arbitrary mul- 
tiplicative constant can be absorbed into C and so (5) can give rise to families of equilibrium 
points. 

If Rank(J(z)) — m — 2n, (4) has only the trivial solution C(z) = constant. If 
Rank(J(z)) = 2n < m and n is constant, then (4) has m — In functionally independent 
solutions (Lie-Darboux theorem). The problem becomes more interesting if there is a sin- 
gularity where Rank(J(z)) changes: in this case we have a singular (hyperfunction) Casimir 
element. For example, consider the one-dimensional system where J = ix (x G R). At 
x = Rank(J) drops to 0, and this point is a singular point of the differential equation 
J(x)d x C = 0. The singular Casimir element is C(x) = Y(x), where Y is the Heaviside step 
function. 

We generalize (1) further to include infinite-dimension systems. Let u G X be the state 
variable, where for now X is some unspecified function space, J{u) be a linear antisymmetric 
operator in X that generally depends on u (for a fixed u, J{u) may be regarded as a linear 
operator X — > X - see Remark 1 below), and H(u) be a functional X — > IBL Introducing an 
appropriate functional derivative (gradient) d u H{u), we consider evolution equations of the 
form 

dtu = J(u)d u H(u) , (6) 

where dtu := du/dt. A Casimir element C(u) (a functional X — > M) is a nontrivial solution 

to 

J(u)d u C(u) = 0. (7) 

We may solve (7) by two steps: 



1. Find the kernel of J~(u), i.e., solve a "linear equation" (cf. Remark 1) 

J{u)v = (8) 

to determine v for a given u, which we write as v(u). 

2. "Integrate" v(u) with respect to u to find a functional C(u) such that v(u) = d u C(u). 

As evident in the above finite-dimension example, step-1 should involve "singular so- 
lutions" if J{u) has singularities. Then, step-2 will be rather nontrivial - for singular 
Casimir elements, we will need to generalize the notion of functional derivative. As men- 
tioned above, the present paper is devoted to such an extension of the notion of Casimir 
elements in infinite-dimensional noncanonical Hamiltonian systems. Specifically, we invoke 
the Euler equation of ideal fluid mechanics as an example, but much carries over to other 
systems since many fluid and plasma systems share the same operator J . In Sec. II, we 
will describe the Hamiltonian form of Euler's equation, which places on a more rigorous 
footing the formal calculations of [6-11]. In Sec. Ill A, we will analyze the kernel of the 
corresponding symplectic operator Jiu) and its singularity. A singular Casimir element 
and its appropriate generalized functional derivative (gradient) will be given in Sec. Ill B. 

The relation between the (generalized) Casimir elements and equilibrium points (station- 
ary ideal flows) will be discussed in Sec. IV. Generalizing (5) to an infinite-dimensional 
space, we may find an extended set of equilibrium points by solving 

d u [H{u) + C{u)\ = 0. (9) 

We note, however, it is still uncertain whether or not every equilibrium point can be obtained 
from Casimir elements in this way. For example, let us consider a simple Hamiltonian 
H(u) = ||m|| 2 /2 (in Appendix A, the Hamiltonian of the Euler equation is given in terms of 
the velocity field u, which here corresponds to the state variable u). Then, d u H{u) = u and 
(6) reads 

d t u = J(u)u. (10) 

The totality of nontrivial equilibrium points is Ker(J'(u)). For v G Ker(j7"(«)) to be char- 
acterized by (9), which now simplifies to u — —d u C(u), we encounter the "integration prob- 
lem," i.e., we have to construct C(u) such that v(u) = d u C(u) for every v(u) G Ker(j7"(w)) - 
this may not be always possible. On the other hand, even for a given C(u), the "nonlinear 



equation" u = —d u C(u) does not necessarily have a solution - in Sec. IV B we will show 
some examples of no-solution equations. While we leave this question (the "integrability" of 
all equilibrium points) open, the present effort shows it is sometimes possible and provides 
a more complete understanding of the stationary states of infinite-dimensional dynamical 
systems. In Sec. V we give some concluding remarks. 

Remark 1 In (6), the operator J (u) must be evaluated at the common u of d u H(u), thus 
J~(u)d u H(u) is a nonlinear operator with respect to u. However, the application of J(u) (or 
J[u)d u ) may be regarded as a linear operator in the sense that J~(u)(av + bw) = aj{u)v + 
bj(u)w (or J(u)d u [aF(u)+bG{u)} = aJ(u)d u F(u) + bJ{u)d u G{u)). Note that J(u)v (for 
v = d u F(u)) is not J{v)v. 

II. HAMILTONIAN FORM OF THE EULER EQUATION 

Investigation of the Hamiltonian form of ideal fluid mechanics has a long history. Its 
essence is contained in Lagrange's original work [12] that described the fluid in terms of the 
"Lagrangian" displacement. Important subsequent contributions are due to Clebsch [13, 14] 
and Kirchhoff [15]. In more recent times, the formalism has been addressed in various ways 
by many authors (e.g. [16-21]). Here we follow the noncanonical Poisson bracket description 
as described in [2, 6]. Analysis of the kernel of J requires careful definitions. For this reason 
we review rigorous results about Euler's equation in Sec. II A, followed by explication of 
various aspects of the Hamiltonian description in Sees. II B II E. This places aspects of 
the noncanonical Poisson bracket formalism of [7, 8, 10, 11] on a more rigorous footing; of 
particular interest, of course, is the symplectic operator J that defines the Poisson bracket. 

A. Vorticity Equation 

Euler's equation of motion for an incompressible inviscid fluid is 

d t u + (u ■ V)u = -Vp (in ft), (11) 

V • u = (in ft) , (12) 

n ■ u = (on <9ft) , (13) 



where Q is a bounded domain in IR n (n = 2 or 3) with a sufficiently smooth (C 2+e -class) 
boundary dQ, n is the unit vector normal to dQ, u is an n-dimensional vector field repre- 
senting the velocity field, and p is a scalar field representing the fluid pressure (or specific 
enthalpy); all fields are real- valued functions of time t and the spatial coordinate x E Q. 
We may rewrite (11) as 

d t u = u x u) — Vp , (14) 

where u) = V x u is the vorticity and p = p + u 2 /2 is the total specific energy. The curl 
derivative of (14) gives the vorticity equation 

d t ijJ = V x (u x w) . (15) 

We prepare basic function spaces pertinent to the mathematical formulation of the Euler 
equation. Let L 2 (Q) be the Hilbert space of Lebesgue-measurable and square-integrable real 
vector functions on Q, which is endowed with the standard inner product (a, b) = f Q dx a ■ b 
and norm ||a| = (a,a) 1//2 . We will also use the standard notation for Sobolev spaces. We 
define 

Ll(n) = {ue L 2 (n) | V • u = 0, n ■ u = 0}, (16) 

where n ■ u denotes the trace of the normal component of u onto the boundary dQ, which is 
a continuous map from {u E L 2 {Vt) \ V • u E L 2 (f2)} to H~ l l 2 (dVl). We have an orthogonal 
decomposition 

L 2 (VL) = Ll(VL)®{V6\ 9eH\Q)}. (17) 

Every u E L 2 (Q) satisfies the conditions (12) and (13), thus we will consider (14) to be an 
evolution equation in the function space L 2 (Q) (cf. Appendix A). 

Hereafter, we assume that the spatial domain has dimension n — 2, and O C I 2 is a 
smoothly bounded and simply connected (genus=0) region [22]. For convenience in formu- 
lating equations, we immerse Q C M 2 in IR 3 by adding a "perpendicular" coordinate z, and 
we write e = V z. 

Lemma 1 Every two-dimensional vector field u satisfying the incompressibility condition 
(12) and the vanishing normal boundary condition (13) can be written as 

u = Vy? x e (18) 



with a single-value function tp such that tp\dn = 0, i.e., 

Ll(n) = {Vy? xe\ipe H%(n)}. (19) 

Proof For the convenience of the reader, we sketch the proof of this frequently-used lemma 
(see e.g. [10, 13, 14, 16, 23, 24] for discussions of various potential representations). Evidently, 
V-(Vy2xe) = V-[Vx (tpe)} = 0, and n- (Vy? x e)\ m = (e x n) ■ Vtp\ m = if tp G Hq(Q). 
Thus, the linear space X = {Vtpxe | tp G Hq(Q)} is contained in L 2 a (Q). And, the orthogonal 
complement of X in L%(Q) contains only the zero vector: Suppose that u G L 2 a {Q) satisfies 

(«, Vtp x e) = Vv? G H*(Q). (20) 

By the generalized Stokes formula, we find (u, Vtp x e) = (e • V x u, tp). Since V x u has 
only the e component, (20) implies V x u = 0. Since u G L%.(Q), we also have V • u = 
and n ■ u\ dn = 0. In a simply connected Q, the only such u is the zero vector. Hence, we 
have (19). | 

Using the representation (18), we may formally calculate 

w = Vxm= (— Atp)e =: cue . 

The vorticity equation (15) simplifies to a single e-component equation [25]: 

d t cu = {cu,)Cuj} (infi), (21) 

where 

{a, 6} = — Va x V6 • e = d y a ■ d x b — d x a ■ d y b, 

and K, is the inverse map of —A with the Dirichlet boundary condition, i.e., /C : u \— > tp 
gives the solution of the Laplace equation 

-Atp = cu (in fi), tp = (on dQ). (22) 

As is well-known, /C : L 2 (Q) — > Hq(Q) D H 2 (Q) is a self- adjoint compact operator. For 
tp G i?Q($l), we define u = —Atp as a member of H" 1 ^), the dual space of Hq(Q) with 
respect to the inner-product of L 2 (Q). The inverse map (weak solution), then, defines 



Theorem 1 We regard the vorticity equation (21) as an evolution equation in H 1 (fi), i.e., 
we consider the weak form: 

(d t u - {u, Ku}, 0) = V0 e H%(n) . (23) 

By the relations <p = Ku>, u = Vy? x e, and uj = use, (23) is equivalent to the Euler equation 
(14) as an evolution equation in L^(Q). 

Proof In the topology of L 2 ff (Q), the Euler equation (14) reads as 

(dtu-ux uj + Vp,v) = VveLl(Q). (24) 

By (17), the left-hand side of (24) reduces to (dtu — u x uj,v). By Lemma 1, we may put 
v = V<j) x e = V x (4>e) with e Hq(Q). Finally, plugging this representation into (24), we 
obtain 

{d t u — u x w, V x (0e)) = (e • V x (d t u — u x w), </>) 

= (^w- {w,y9},0). 

Hence, (24) is equivalent to (23). | 

B. Hamiltonian 

Now we consider the Hamiltonian form of the vorticity equation (21) - to be precise, its 
"weak form" (23) (cf. Appendix A which treats the Euler equation of (14)). 

First we note that the natural choice for the Hamiltonian is H = ||w|| 2 /2, the "energy" 
of the flow u. Using u = V(p x e, we may rewrite H = || Vy?|| 2 /2 = (<^>, —Atp)/2. Selecting 
the vorticity u> as the state variable, we define (by relating ip = Kuu) 

H(u) = - J dx (Km>) ■ uj , (25) 



2jq 

which is a continuous functional on iJ _1 (fi). This is equivalent to the square of the norm 
of H~ l (Q), i.e., the negative norm induced by Hq(Q). 

C. Gradient in Hilbert Space 

Next, we consider the gradient of a functional in function space. Let $(w) be a functional 
defined on a Hilbert space X. A small perturbation eu G X (|e| ^C 1, ||u|| = 1) will 



induce a variation 8Q(u;u) = $(w + eu) — <&(u). If there exists a g G X* = X such that 
5$(w; -u) = e(g, u) + 0(e 2 ) for every u, then we define (9 u $(m) = g, and call it the gradient of 
$(«). Evidently, the variation |5$(-u; u)\ is maximized, at each u, by u = d u <&(u)/\\d u <&(u)\\. 
The notion of gradient will be extended for a class of "rugged" functionals, which will be 
used to define singular Casimir elements in Sec. Ill B. As for the Hamiltonian, however, we 
may assume it to be a smooth functional. The pertinent Hilbert space is L 2 (Q), on which 
the Hamiltonian H(u) is differentiable; using the self-adjointness of /C, we obtain 

d^Hiuj) = Kuj. (26) 

Note that the gradient <9 w if(co>) may be evaluated for every uo G H~ l (Q) with the value in 

D. Noncanonical Symplectic Operator 

Finally, we describe the noncanonical symplectic operator J(uj) of [7, 8, 10, 11]. Formally, 
we have 

J{u)i) = \{d y uj)d x - (d x co)d y ] i) = {u, ^} , (27) 

which indicates u must be a "differentiable" function. However, we will need to reduce this 
regularity requirement on u. Thus, we turn to the weak formulation that is amenable to the 
interpretation of the evolution in if _1 (f2) (see Theorem 1). Formally, we may calculate 

(J(u>)1>, 0) = ({co, ^}, 0) = (to, {i>, 0}) , (28) 

with the right-hand-side finite (well-defined) for u G C(Q) and ijj,(f) G Hq(Q). In fact, 

\(u,{ip,(j>})\ < llwllsup dx\{i>,(f)}\ 

Jn 

< |M| SU p [dx\Vi/>\\V(i>\ 

< |M| sup ||W||||V0||, 

where IMIsup = su P x eQ.\ U! ( x )\- Hence, we may consider the right-hand-side of (28) to be 
a bounded linear functional of G Hq(Q), with uj and ip acting as two parameters. We 
denote this by (ou, {-0,0}) =: F(w,-0;0). By this functional, we "define" J{oj)ip on the 
left-hand-side of (28) as a member of Hq(Q)* = 7J _1 (f2), i.e., we put 

(J(u)rJ>,<j>):=F(u,rl>;<j>) V0 G H*(Q). 



For a given u G C(Q), we may consider that J~(oj) is a bounded linear map operating 
on V, i.e., J(u) : #o(fi) ->■ # _1 (Q). Evidently, {J{u)^,<f>) = -($ , J(u)<j>) , i.e., J"(w) is 
antisymmetric. 

E. Hamiltonian Form of Vorticity Equation 

Combining the above definitions of the Hamiltonian H(u), the gradient d u , and the 
noncanonical symplectic operator J7"(w), we can write the vorticity equation (21) in the 
form 

&tu = J{u)d u H{u). (29) 

As remarked in Theorem 1, (29) is an evolution equation in if _1 (f2) (cf. Appendix A for the 
u formulation). 

For every fixed cu G C(Q), J{u) may be regarded as a bounded linear map of Hq(Q) — > 
iy _1 (fi), where the bound changes as a function of u. And d u H(u) is a bounded linear 
map of i/ _1 (r2) —?• Hq(Q). Hence, each element composing the right-hand side (generator) 
of the evolution equation (29) is separately regular. However, their nonlinear combination 
can create a problem: As noted in Remark 1, we must evaluate the operator J{u) at the 
common u> of d w H(u). While d u H{u) can be evaluated for every ui G H~ l (Q) with its 
range = Hq(Q) = domain of J~(uj), if J{oj) is defined; however, we can define the operator 
J{l>j) only for uo G C(Q). The difficulty of this nonlinear system is now delineated by the 
singular behavior of the symplectic operator J{oo) as a function of u - if the orbit u(t) (in 
the function space ff _1 (f2)) runs away so as to increase ||o;|| SU p, the evolution equation (29) 
will breakdown. 

To match the combination of J{uj) and d u H{tjj), the domain of the total generator 
J (10)8^11 (u) must be restricted in C(Q). Fortunately, this domain is not too small; the 
regular (classical) solutions for an appropriate initial condition lives in this domain, i.e., if 
a sufficiently smooth initial condition is given, the orbit stays in the region where ||w|| S up is 
bounded [26]. 
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III. CASIMIR ELEMENTS 

A. The Kernel of J(u) 

We begin with a general representation of the kernel of the noncanonical symplectic 
operator J~(uj), which will be a subset of its domain Hq(Q) (see Sec. II D). 

Proposition 1 For a given u> G C(Q), ip G Hq(Q) is an element of Ker (J(u})), iff 

ooVip = V6 36 G H\Q), (30) 

which implies that 

Kei(J(u)) = {*pe H^{Q) | uV^ G L^ft)" 1 }. (31) 

Proof By the definition (27), if) G Ker(j7"(cu)) implies {u, if)} = in the topology of il _1 (f2), 
i.e., 

({u, if)}, 0) = -(uVip, V0 x e) = V0 G H*(n). (32) 

By Lemma 1, (32) implies that 

(uVi«) = V?; G L^(fi). (33) 

Remembering (17), we obtain (30) and (31). | 

To construct a Casimir element from if) G Ker(j7"(o;)), we will need a more "explicit" 
relation between u and ■0. We will show how such a relation is available for a sufficiently 
regular u. 

Let us start by assuming uj G C l (VL). Then, we may evaluate J{lS)i\) as {oo,if>} = 
-Vw xVf eGl 2 (0). Therefore, ^ G H%(Q) belongs to Ker (J (to)), iff 

{co,iP} = [GL 2 ((])]. (34) 

Equation (34) implies that two vectors Vw G C($l) and V^ G L 2 (Q) must align almost 
everywhere in Q, excepting any "region" in which one of them is zero. Such a relationship 
between u and if) can be represented, by invoking a certain scalar ((x,y), as 

w = /(0, V> = <7(0- (35) 
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FIG. 1: Depiction of Q (uj), the largest subset of Q bounded by an uj level set. Note, the complement 
of Q {ui) contains level sets that intersect dQ.. Also depicted is the plateau set Q p C Q, a region 
where w=constant. 

The simplest solution is given by tp = g{oS) (i.e., / = identity). In later discussion, we shall 
invoke a nontrivial / to represent a wider class of solutions. 

We note that the condition ip G Hq(Q) implies the boundary condition ip\gn = 0. For 
(34) to be compatible with this boundary condition, = n • (Vcu x V^) = (n x Vw) • Vip on 
dQ. Hence, to find ip such that Vip ^ on dQ, we need w|an = constant. Otherwise (i.e., 
if u ytz constant on some dQ' C dQ), ip must vanish in a neighborhood of dQ' . To formulate 
the solvability (compatibility) condition, we denote by Q (u) the largest region in Q (not 
necessarily a connected set) which is bounded by a level set (contour) of u. See Fig. 1. 
If w|ao 7^ constant, Q (u>) is smaller than Q, and then, every level set of u in Q \ Q (u) 
intersects dQ. Hence, supp(?/>) := closure of {x G Q \ if)(x) ^ 0} C Q (uj). We shall assume 
that Q {oj) 7^ for the existence of nontrivial ip G Ker(j7"(a;)). 

Now we make a moderate generalization about the regularity: Suppose that uj is Lipschitz 
continuous, i.e., uj G C 0,1 (Q). Then, uj has a classical gradient Vw almost everywhere in Q 
(see e.g. [27]), and |Vw| is bounded. Note, uj may fail to have a classical Vuj on a measure- 
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zero subset Q s C f2, but for x G fi s we may define a set-valued generalized gradient (see 
[28]). With a Lipschitz continuous function g : 1 -> 1 we can solve (30) by 

ij> = g(u>), 9 = &(u), (36) 

where $(£) = <?(£)£ — fd£g(£). To meet the boundary condition V'lan = 0, g must satisfy 

£(w(a;)) = Vx e 9Q. (37) 

However, the solution (36) omits a different type of solution that emerges with a singu- 
larity of J~(oj): If w has a "plateau," i.e., cu = Uq (constant) in a finite region Q p C Q (see 
Fig. 1), the operator J(oj) trivializes as J{oj) = {u , •} = in Q p (i.e., the "rank" drops to 
zero; recall the example of Sec. I), and within Q p we can solve (30) by an arbitrary if) with 
6 = UQip. Notice that the solution (36) restricts if) = g(uo) = constant in Q p . To remove 
this degeneracy, we abandon the continuity of g and, for simplicity, assume that u has only 
a single plateau. First we invoke the reversed form [cf. (35)]: 

u = f(if)) , (38) 

where we assume that / is a Lipschitz continuous monotonic function. Denoting $(77) = 
jdi] f(i]), we may write uVif) = Vi?(V0, where the gradients on both sides evaluate classically 
almost everywhere in Q if if> is Lipschitz continuous. If the function f(ip) is flat on some 
interval, f{ip) = uq = constant for ip_ < ip < ip + , a plateau appears in the distribution of u. 
See Figs. 1 and 2(a). Since the present mission is to find ip for a given u, we transform (38) 
back to (36) with the definition g = f~ l . (In Sec. IV, however, we will seek an equilibrium u 
that is characterized by a Casimir element, and then, the form (38) will be invoked again). 
A plateau in the graph of / will, then, appear as a "jump" in the graph of g. See Fig. 2(b). 

We now allow the function g(oo) to have a jump at uq = oo\q p . Formally we write 
g(uj) = ^l(w) + ctY{u — ujo), with gL,{oj) a Lipschitz continuous function, Y(u> — uj ) the step 
function, and a a constant determining the width of the jump. We connect the graph of the 
step function by filling the gap between lim w -|- wo g(u>) = ■0- and lim^ g(u>) = ip+ = i ) - +«• 
See Fig. 2(b). Since g(u) is multi- valued at u = Uq, ip(x) = g(u(x)) is arbitrary in the range 
of [if)-,ip + ], i.e. in Q p . Choosing a sufficiently smooth ip in Q p , we may assume ip G H 1 ^). 

Summarizing the above discussions (and making an obvious generalization), we have 
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FIG. 2: Depiction of the dual relationships between ip and u. (a) The relation uj = f(ip) with 
the plateau singularity, (b) The singular (discontinuous) function ip = g(oj) illustrating the kernel 
element that stems from a plateau of u. 

Corollary 1 Suppose that u G C 0,1 (Q) andQ (u) ^ 0. ThenKer{J{u:)) contains nontrivial 
elements, and a part of them can be represented as 



V> = </(w), 



(39) 



where g(£) is an arbitrary function that satisfies the boundary condition (37) and such that 

V 

9(0 = 9L(S) + J2a t Y{S-u> t ), (40) 



e=i 



where ug denotes the height of a plateau ofu, Y(£) is the 'filled" step function, and a^ is a 
constant. 

Remark 2 Clearly the form (40) of g(u) is rather restrictive: 

(i) In the plateau region, (30) has a wider class of solutions. In fact, ip may be an 
arbitrary H^^-class function whose range may exceed the interval [ip-,ip + ]. In this case, the 
graph of g(u>) has a "thorn" at u , and we may not integrate such a function to define a 
Casimir element G(u). (See Sec. IIIB.) 

(ii) In (40), we restrict the continuous part gi{0 to be Lipschitz continuous, by which 
i\) = g{uS) (00 G C 0,l (Vt)) is assured of Lipschitz continuity (thus, ip G H l {VL)). However, 
this condition may be weakened, depending on the specific u, i.e., it is only required that 
g'i^VueH^Q). 
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B. Construction of Casimir Elements 

Our next mission is to "integrate" the kernel element ip G Ker(j7"(cu)) as a function of 
u, and define a Casimir element C(ui), i.e., to find a functional C(u) such that <9 w C(w) G 
Ker(j7"(w)). To this end, the parameterized ip of (39) will be used, where the function 
g(u) may have singularities as described in Corollary 1. The central issue of this section, 
then, will be to consider an appropriate "generalized functional derivative" by which we 
can define "singular Casimir elements" pertinent to the singularities of the noncanonical 
symplectic operator J {us). 

Let us start by considering a regular Casimir element generated by g(£) G C(R): 

C G (u>)= fdxG(u), (41) 

Jn 

where G(£) := fd£g(£). The gradient of this functional can be readily calculated with the 

definition of Sec. II C: Perturbing u by eu results in 

5Cg(u, oo) — e / dx g{uo)uj + 0(e 2 ) . 
Jn 

Hence, we obtain <9 w Cg(w) = g(u>), proving that Ca(oo) of (41) is the Casimir element 

corresponding to g(uj) G Ker(j7"(o;)). 

Now we construct a singular Casimir element corresponding to a general g(u) that may 
have "jumps" at the singularity of J~{oo), i.e., the plateaus of u. The formal primitive 
function G(£) of such a g(£,) has "kinks" where the classical differential does not apply 
— this problem leads to the requirement of an appropriately generalized gradient of the 
functional Cg{<^) = J n dxG(ou) generated by a kinked G(£). 

Here we invoke the Clarke gradient [28], which is a generalized gradient for Lipschitz- 
continuous functions or functionals. Specifically, if F : R — > R, then the Clarke gradient of 
F at x, denoted by d x F(x), is defined to be the convex hull of the set of limit points of the 
form 

lim d x F(x + 8j) with lim Sj = 0. (42) 

Evidently, d x F(x) is equivalent to the classical gradient, d x F(x), if F(x) is continuously 
differentiable in the neighborhood of x. Also, it is evident that a "kink" in F yields d x F(x) 
with a graph that has a "jump" with the gap filled as depicted in Fig. 2(b). When F{u) 
is a convex functional on a Hilbert space X, i.e. F : X — > R, then d u F(u) is equal to the 
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sub- differential: 

d u F:u^{g\ F(u + 5)- F(u) > (g, u), V5 e X}, (43) 

which gives the maximally monotone (i.e., the gap-filled) function [29]. For the purpose of 
Sec. IV, a monotonic (?(£) = c^G(£) will be sufficient. 

From the above, the following conclusion is readily deducible: 

Proposition 2 Suppose that uj e C°' l (Q) andQ (u) ^ 0. By g(£) satisfying (37) and (40), 
we define G(£) such that g(£) = d^G(^). Then, Cg(u>) = J n dxG(u) is a generalized Casimir 
element, i.e., d^Ccijjj) G Ker[J{ui)). 

IV. EXTREMAL POINTS 

A. Extremal Points and Casimir Elements 

The Casimir element Cg{oo) naturally extends the set of extremal equilibrium points, 
the set of solutions of the dynamical system (29) that are both equilibrium solutions and 
energy- Casimir extremal points. This is done by including extremal points satisfying [cf. 

(9)] 

B u [H(u) + C G (u)]3 0, (44) 

which explicitly gives 

Kmj E -g{oj). (45) 

Here we assume that <?(£) = d%G(£) is a monotonic function, i.e., G(£) is convex. Then, 
we may define a single-valued continuous function / = (— g)^ 1 and rewrite (45), denoting 
ip = 1Cu>, as 

-Atp = f(<p). (46) 

If /(0) 7^ 0, (46) will determine a nontrivial (tp ^ 0) equilibrium extremal point. Here we 
prove the following existence theorem: 

Theorem 2 There exists a finite positive number M, determined only by the geometry of 
Q, such that if 

L ;= inf ^M < M (47) 

V £R \r]\ 

then (46) has a solution (p G Hq(Q). 
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Proof We show the existence of the solution by Schauder's fixed-point theorem (see for 
example [30] and [31, p. 43]). First we rewrite (46) as 

u = F(u) ■= f(Kw), (48) 

where /C is a compact operator on L 2 (Q) (cf. (22)). Since / G C(R), J 7 is a continuous 
compact map on L 2 (Q). We consider a closed convex subset 

W K = {u\ \\lo\\ < A} Ci 2 (0), 

where ||w|| is the norm of L 2 (Q) and the parameter A will be determined later. We show 
that the compact map J 7 has a fixed point in W\. By Poincare's inequality, we have, for 

y? G H^n), 

IMI<ci||v^||, 

where c\ is a positive number that is determined by the geometry of Q. For u = Vy? x e G 
LKtyDH^tt), we have 

|| ix || < c 2 || V x w||, 

where C2 is also a positive number determined by the geometry of Q. Here ||u|| = || S/ip\\ and 
|| V x n|| = ||A(/?||. Hence, denoting c p = c<iC\ + c-i + 1, we have 

IMI^:=IMI + HVy?|| + ||Av?|| <c p ||A^||. 

By Sobolev's inequality, we have 

sup \tf\ < C S ||<^|| H 2, 

n 
where c s is again a positive number determined by the geometry of Q. Summarizing these 
estimates, we have, for u> G W\, 

sup | /Co; | < c s Cp||a;|| < c s c p A. 
n 

For an arbitrary positive number S, there is a finite number rj* such that |/(±?7*)| < (L + 

5) 1 1]* | . If we choose A = \r)*\/(c s c p ), then by the monotonicity of f(i]), we find 

sup |/(/Cw)| < (L + <J)|t7*| = (L + 5)c s c p A, (49) 

a 

and, upon denoting |f2| := f n dx, we obtain 

||/(/Cw)|| < \Vt\ 1/2 {L + 5)c s c p A= (L + S)A/M, 
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where 

M:= wh~/ (50) 

If (L + 5) < M, we estimate ||/(/Cw)|| < A, and thus J 7 maps W\ into W\. Therefore, T 
has a fixed-point in W\. If W\ contains a fixed point, then an arbitrary W\> with A' > A 
contains a fixed-point. Since S is arbitrary, L < M is the sufficient condition for the existence 
of a fixed point. § 

Corollary 2 7/L = 1/(77*) |/|t/*| > for some finite 77* , £/ien we mm/ assume 5 = in (49), 
and the solvability condition (47) can be modified to be L < M. 

The bound M defined by (50) is related to the eigenvalue of —A = /C _1 . On the other 
hand, the constant L defined by (47) is a property of the Casimir element. As noted above, 
because of the homogeneity of (4), there is at least a one-parameter family of Casimir 
elements and the choice of a multiplicative constant reciprocally scales L. This provides an 
arbitrary parameter multiplying the Casimir element in the determining equation (44), and 
such a parameter may be considered to be an "eigenvalue" characterizing the stationary 
point. The no-solution example of the next section will reveal the "nonlinear property" of 
this eigenvalue problem. 

B. No- Solution Example 

Here we present a class of examples that violate the solvability condition proposed in 
Theorem 2. This demonstrates that finding Casimir elements that generalize the extremal 
equation (46), does not automatically extend the set of extremal equilibria, since the result- 
ing elliptic equation may not have a solution. This is demonstrated by the following two 
propositions. 

Proposition 3 Let f(tp) have the form 

f(<p)=\i<p + &(<p), (51) 

where Ai is the first eigenvalue of —A = /C _1 and $ > 0. Then, (46) does not have a 
solution. 
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Proof We denote by <pi the eigenfunction corresponding to the eigenvalue Ai. Upon taking 
the inner product of (46) with 0i, we obtain 

(-Atp, 00 = (A lV 9, 00 + (#(tp), 0i). (52) 

The left-hand-side of (52) satisfies (— A</?, 0i) = (ip, — A0i) = (Ai</?, 0i), which cancels the 
first term on the right-hand-side, leaving (^(<^),0 1 ) = 0. However, this is a contradiction, 
because 0i > (or < if otherwise normalized) on Q (cf. [32]) and, by the assumption, 
*(¥>) > o. I 

The above result is easily generalized as follows: 

Proposition 4 Suppose that 

f(<p) = \ s <p + #(<p), (53) 

where Xj is any eigenvalue of —A and < A < $((p) < B, for some A and B. Let <pj be the 
eigenfunction corresponding to Xj. Since <j)j may not have a definite sign in Q, we divide Q 
into Q + = {i G H <Pj(x) > 0} and Q_ = {i G O <f>j(x) < 0}, and define 



P+= dxfa, P =- dx fa (both > 0) 
If P + ^ P_, (46) does not have a solution for some ranges of A and B. 
Proof As shown in the proof of Proposition 4, a solution of (46) must satisfy 



(ti(tp),<j>j)= / dx^(if)(j) j + / dx$((p)(j) j = 



By assumption, we have 

A(j)j < -&((p)(f)j < B(f)j in Q + , 
A(f>j > f?(<^)0j > B(j)j in fi_. 

Hence, the inequalities 

AP + - BP_ < < BP + - AP_ 

must hold, i.e., 

A/B < P./P+, (54) 

B/A > P-/P+. (55) 

However, if P_/P + < 1, and A and B are such that 1 > A/B > P_/P + , then there is a 
contradiction with (54). Similarly, if P_/P + > 1, and A and B are such that 1 < B/A < 
P-/P + , then there is a contradiction with (55). | 
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V. CONCLUDING REMARKS 

After establishing some mathematical facts about the Hamiltonian form of the Euler 
equation of two-dimensional incompressible inviscid flow, we studied the center of the Lie- 
Poisson algebra, i.e., the kernel of the noncanonical symplectic operator J{oo). Casimir 
elements C{u) were obtained by "integrating" the "differential equation" 

J{w)d u C(u) = . 

For finite-dimensional systems with phase space coordinate z, this amounts to an analysis 
of V := J{z)d Z) a linear partial differential operator, and nontriviality can arise from a 
singularity of V, whence an inherent structure emerges. Recall the simple example given in 
Sec. I where V = ixd x was seen to generate the hyperfunction Casimir C(x) = Y(x). For 
finite-dimensional systems the theory naturally finds its way to algebraic analysis: in the 
language of D-module theory, Casimir elements constitute Ker(P) = Horrid (Coker('P), F), 
where T> is the ring of partial differential operators and F is the function space on which V 
operates, and Coker('P) = T> jT>V is the D-module corresponding to the equation VCiuS) = 0. 
However, in the present study u is a member of an infinite-dimensional Hilbert space, thus 
V may be regarded as an infinite-dimensional generalization of linear partial differential 
operators. From the singularity of such an infinite-dimensional (or functional) differential 
operator V = J~(oo)d w , we unearthed singular Casimir elements, and to justify the operation 
of V on singular elements, we invoked a generalized functional derivative (Clarke differential 
or sub-differential) that we denoted by d w . 

For infinite-dimensional systems, we cannot "count" the dimensions of Ker(P) and 
Ker(J'). It is, however, evident that dim-Ker('P) < dim-Ker(j7"), if J has singularities, 
i.e., singularities create "nonintegrable" elements of Ker(J'). As shown in Corollary 1, a 
plateau in u> causes a singularity of J{ui) and generates new elements of Ker( l 7(w)), which 
can be integrated to produce singular Casimir elements (Proposition 2). However, as noted 
in Remark 2(i), more general elements of Ker(j7"(cu)) that are not integrable may stem from 
a plateau singularity. Moreover, we had to assume Lipschitz continuity for uj to obtain an 
explicit relation between ip e Ker(j7"(o;)) and u - otherwise we could not integrate ip with 
respect to u to construct a Casimir element. In the general definition of J~(oj), however, u 
may be nondifferentiable (we assumed only continuity), and then, a general ip G Ker(j7"(w)) 
may not have an integrable relation to u (see Proposition 1). 
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In Sec. IV we solved the equation 

B u [H(u) + C(u)] 3 

for uo, where the solution u gave an equilibrium point of the dynamics induced by a Hamilto- 
nian H(u). A singular (kinked) Casimir yielded a multivalued (set- valued) gradient d u C(u), 
encompassing an infinite-dimensional solution stemming from the plateau singularity. This 
arose because in the plateau of u, ip G Ker(j7"(w)) is freed from u and may distribute arbi- 
trarily. The component 5'l(w) of the Casimir C(u) of (40) represents explicitly the regular 
"dimensions" of Ker(j7"(o;)). In contrast, the undetermined dimensions pertinent to the 
singularity ui = ujq are "implicitly" included in the step-function component of (40), or in 
the kink of C(u). However, for a given Hamiltonian H(u), i.e. a given dynamics, a specific 
relation between ip G d u C{uj) and uj emerges. 

Theorem 2 of Sec. IV A and the nonexistence examples of Sec. IV B may not be new 
results in the theory of elliptic partial differential equations, but they do help delineate 
the relationship between Hamiltonians and Casimir elements, viz. that Casimir elements 
alone do not determine the extent of the set of equilibrium points. In the present paper, 
we did not discuss the bifurcation of the equilibrium points; the reader is referred to [31] 
for a presentation of the actual state of the studies of semilinear elliptic problems and 
several techniques in nonlinear analysis (see also [30] and [33]). We also note that we have 
excluded nonmonotonic g(oo) that will make f((fi) multivalued or, more generally, equations 
like $(A(/?, (p) = 0; cf. (35). For fully nonlinear elliptic partial differential equations, the 
reader is referred to [30, 34-37]. 
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Appendix A. On the symplectic operator in terms of the velocity field u G L 2 (Q,) 

Here we formulate the Euler equation, for both n = 2 and 3, as an evolution equation in 
L 2 (Q) (see Sec. II A), and discuss the symplectic operator J[u) in this space. This differs 
from the formulation of Sec. II in that the state variable here is the velocity field u G L 2 (Q) 
instead of the vorticity u. As noted in Sec. II A, L 2 (Q) is a closed subspace of L 2 (Q), and we 
have the orthogonal decomposition (17). We denote by P a the orthogonal projection onto 
L 2 (Q). Upon applying P a to the both sides of (14), we obtain 

d t u = —P a iy x u) x u, (56) 

which is interpreted as the evolution equation in L^(f2), where the incompressibility condi- 
tion (12) and the boundary condition (13) are implied by u G L 2 (Q). 
For u G L 2 (Q), the Hamiltonian is, of course, the kinetic energy 

H(u)= 1 -\\u\\ 2 . (57) 

Fixing a sufficiently smooth u acting as a parameter, viz. V x u G C(Q), we define for 
v G L 2 (Q) the following the noncanonical symplectic operator: 

J(u)v = -P CT (V x u) x v . (58) 

As a linear operator (recall Remark 1 of Sec. I), J{u) consists of the vector multiplication 
by (V x u) followed by projection with P a . Evidently, J(u) is antisymmetric: 

(J(u)v, v') = -(v, J{u)v') Vi>, v' G L 2 ff (n). 

In fact, for every fixed smooth u, ijiu) is a self-adjoint bounded operator in L 2 (Q). 

With the H{u) of (57), the J{u) of (58), and the gradient d u (see Sec. II C), we may 
write (56) as 

d t u = J(u)d u H(u). (59) 

Assuming n = 2, which we do henceforth, by Lemma 1 we may put u = Vy? x e and 
u = — A(/9 with tp G Hq(Q). Fixing u G C(Q) as a parameter, and putting v = Vip x e with 
ip G Hq(Q), we may write 

J(u)v = -P a [toe x (W x e)] = -P a [wV^] . 
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By (17), v E Ker(J(u)) iff 

ooVip = V6 BOeH 1 ^), (60) 

which is equivalent to (30). Arguing just like in Sec. Ill A, we find solutions of (60) of the 
form 

ij} = g{u). (61) 

The Casimir element constructed from (61) is 

C G {u)= I dxG{e-V x«)= fdxG(u). (62) 

Jn Jn 

Perturbing u by eu and restricting n ■ u = on dQ yields Su = ee ■ V x u, and 



5C G (u; u) = e dx G'(u)e ■ V x u + 0(e 2 
Jn 



= e / dxVG'(u) xe-u + 0(e 2 ). 
Jn 

Hence, upon denoting g(ui) = G'(u) we obtain 

d u C G (u) = Vg(uj)xe. (63) 

By (61), it is evident that J{u)d u Cc(u) = — P a [uiV g{oj)] = 0, confirming that d u Cc(u) G 
Ker{J(u)). 

Since here, H(u) = ||u|| 2 /2, we may find an equilibrium point by solving u = —d u CG{u) 
(c.f. Sec. I, Eq. (9)). 
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